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We present a multimode theory of non-Gaussian operation induced by an imperfect on/ofï-type 
photon detector on a splitted beam from a wideband squeezed light. The events are defined for finite 
time duration T in the time domain. The non-Gaussian output state is measured by the homodyne 
detector with finite bandwidh B. Under this time- and band-limitation to the quantm states, we 
develop a formalism to evaluate the frequency mode matching between the on/off trigger channel 
and the conditional signal beam in the homodyne channel. Our formalism is applied to the CW and 
pulsed schemes. We explicitly calculate the Wigner function of the conditional non-Gaussian output 
state in a realistic situation. Good mode matching is achieved for BT < 1, where the discreteness 
of modes becomes prominant, and only a few modes become dominant both in the on/off and the 
homodyne channels. If the trigger beam is projected nearly onto the single photon state in the 
most dominant mode in this regime, the most striking non-classical effect will be observed in the 
homodyne statistics. The increase of BT and the dark counts degrades the non-classical effect. In 
the CW scheme, one will be able to attain a stringent mode matching in the regime BT < 0.1. 
Actually the bandwidth can be set typically to B « 10 MHz by using an appropriate set of filters, 
while the duration T can be set in 10 ns order by electrical gating. The spatial mode matching can be 
fulfiUed by careful cavity locking. The temporal mode matching is not a problem in that time scale. 
So the CW scheme will provide a good test-bed for non-Gaussian operations. In the pulsed scheme, 
the duration T is automatically set by the làser pulse width, which is typically of ps order. The 
band limitation comes from the LO bandwidth « 1 THz, corresponding to the Fourier-transform 
limit BT « 1. This still satisfies the condition for observing the non-classical effect provided that 
the mode matching in terms of the other degrees of freedom is perfect. In practice, however, it may 
be more challenging to realize a high quality spatiotemporal mode matching in the pulsed scheme 
than the CW setting. 
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I. INTRODUCTION 

Homodyning and photon counting are Standard tech- 
niques to measure quantum states of optical fields. The 
former accesses the continuous nature of optical fields, 
while the latter does the discrete nature of them. The ho- 
modyning technique can now have a quantum efliciency 
exceeding 99% at certain wavelengths, and has been suc- 
cessfully applied to various tasks in quantum òptics and 
quantum Information science. In particular, it provides a 
powerful tool to completely characterize a quantum state 
by reconstructing the Wigner function of the density op- 
erator, referred to as the quantum state tomography. 

Photon counting technique, on the other hand, still 
fails in attaining near-unit quantum efficiency and single 
photon resolution. These properties, however, are not al- 
ways necessary. In fact the on/ofï type photon counting 
technique with less than unit quantum efficiency can be 
useful for certain applications. One of such applications 
is the conditional photon substraction from the squeezed 



light. That is, a small fraction of the squeezed light 
is beamsplitted as trigger photons, and guided into the 
on/ofï detector. By making the amount of fraction of the 
trigger beam small enough, the detection click effectively 
realizes the single photon subtraction from the main 
signal beam. This method allows one to conditionally 
generate the Schròdinger-cat-like state Q, and also to 
conditionally increase the entanglement of the bipartite 
squeezed beams |^ |^ Q . More importantly these condi- 
tional operations are non-Gaussian operations that can 
be implemented with current technologies. The impor- 
tance of non-Gaussian operations on continuous variables 
(CVs) is highlighted by the foUowing no-go statements 
concerning the Gaussian operations. Firstly, quantum 
speed up is impossible for harmònic oscillators by Gaus- 
sian operations with Gaussian inputs 0. Secondly, the 
distillation of Gaussian entanglement from two Gaussian 
entangled states is impossible onlyby Gaussian LOCC 
based on homodyne detection [aQjlS· Therefore, im- 
plementation of non-Gaussian operations is crucial to ex- 
tract ultimate potential of photonic quantum Information 
processing. 
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The first experimental demonstration of the 
measurement-induced non-Gaussian operation was 
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done by Wenger et al. in the ultra-short pulsed regime 
01 . They could observe the phase sensitive non-Gaussian 
distributions, that is, a small dip around the origin of 
the phase space in an unisotropic Wigner function 
distribution. In order to observe the negative dip in 
the Wigner function distribution, which is a strong 
indication of the non-classicahty of the state, more 
careful mode matching considerations should be taken. 

Ideally the main signal beam and the trigger beam 
must be prepared in the same mode. That is, trigger 
photons detected by the photon counter must be in a 
mode obseved by the homodyne detector with respect to 
the spatial, temporal, frequency and polarization modes. 
Otherwise the homodyne detector will see the modes that 
are not quantum correlated with the trigger photons, de- 
grading the non-Gauusian operations significantly. 

The mode matching problem in this kind of condi- 
tional operation was first studied by Ou fï^ . He studied 
the mode matching in the scheme of homodyne measure- 
ment of a conditionally prepared single photon state by 
measurements on a biphoton state produced in the para- 
metric down-conversion. This scheme was originally pro- 
posed by Yurke and Stoler [ïï|, which was analyzed in 
a single mode assuming the perfect mode matching. Ou 
studied it in more practical situations, and showed that 
one has to use a narrow spectral filter in the trigger chan- 
nel in order to match the conditionally prepared single 
photon mode to the one of the LO. 

Grosshans and Grangier extended the analysis, and 
provided a useful formula to renormalize the time and 
frequency overlap between the signal and trigger wave 
packets into an effective quantum efficiency of the homo- 
dyne detector (IJ]. They studied the two cases: (i) con- 
tinuous experiment where the pump and the LO beams 
are initially continuous wave (CW) with small enough 
linewidths, and (ii) pulsed experiment where pump pulses 
short enough are used so that their linewidth becomes 
greater than the frequency filter bandwidth for the trig- 
ger photons. They concluded that the mode match- 
ing condition is much more easily fulfilled in the pulsed 
regime. 

The mode matching in the pulsed regime was further 
studied by Aichele et al. including more general models 
for the spatial and spectral filters in the trigger channel 
[ï^ . They performed an explícit calculation of the degree 
of mode matching in terms of both frequency-momentum 
and time-space representations, and showed ways to at- 
tain the optimal mode matching by using narrowband 
filters in the trigger channel. Viciani et al. further in- 
vestigated temporal and spectral properties of entangled 
photon pair for various filter functions |l4j. 

In these works mm m the parametrlc down 
conversion process is treated within the first order per- 
turbation theory, and it is assumed that the idler photon 
is detected by an ideal photon counter. The conditional 
signal state is then measured by an imperfect homodyne 
detector with the effective quantum efficiency. In the ex- 
periment by Wenger et al., however, the squeezed state 



from a parametric amplifier was used, which is beyond 
the first order perturbation theory of parametric process. 
The state includes higher photon numbers. When the 
trigger beam from such a state is detected by a realistic 
photon counter which can detect the arrival of photons 
but usually fails in precise discrimination of the photon 
number, the conditional output state results in a mixed 
state. Realistic photon counters also suffer from dark 
counts, which cause fake triggers, and degrade the qual- 
ity of the conditional state. 

These imperfections were considered in the context of 
observing the negativity of Wigner function distribution 
of a photon-subtracted squeezed state in The on/off 
detector characteristics including the dark counts was 
modeld by the modal purity factor, and the mode mis- 
match by the inefficiency of homodyne dtector. Unfortu- 
nately, however, theories at such a level of phenomeno- 
logical parameterization cannot provide practically use- 
ful design guidelines. In fact, experimentalists want to 
know, for example, an appropriate time duration of pho- 
ton counting and homodyne detection in order to observe 
the negativity of Wigner function distribution, depend- 
ing on the spectral characteristics of the squeezed state 
sourcc. 

In this paper we develop a multimode theory that can 
explicitly calculate the photon-subtracted squeezed state, 
and can clarify the conditions to observe the negativ- 
ity of Wigner function distribution in realistic situations, 
taking practical imperfections and multimode aspects of 
the quantum states into account. Particular emphasis is 
made on the frequency mode matching under the time- 
and band-limitation to the optical fields. Analysis along 
this linc is essential for designing the practical setup of 
the non-Gaussian operations and analyzing the measure- 
ment results. Nevertheless explict analyses have never 
been given so far to our knowledge. In fact, this is 
the most non-trivial aspect in the mode matching issue 
This problem is also related with the simultane- 
ous control of discrete and continuous natures of quan- 
tum optical fields. Discrete aspects such as photon count- 
ing are often seen in the time domain, while continuous 
ones have been exploited most so far in the frequency do- 
main. The optimal control of both aspects relates deeply 
to the optimal mode matching of the time- and band- 
limited quantum states. The necessary theoretical basis 
was already given by Zhu and Caves [Ï3I . We extend this 
theory into the scheme of conditional preparation of non- 
Gaussian state by photon subtraction from the squeezed 
state. 

In what foUows, the mode matching in terms of the 
other degrees of freedom, i.e. the spatial and polariza- 
tion degrees of freedom, is neglected for simplicity. It is 
actually made almost perfect in the CW scheme based 
on cavity optical parametric oscillator systems by careful 
cavity lockings. We start with a brief review of bàsic no- 
tions and mathematical tools on time- and band-limited 
signals in Section^l Using this basis, we then develop 
a formulation to describe and analyze the mode that the 
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homodyne detector actually observes in Section lïïïl In 
Section lIVl we model the on/ofF detector including prac- 
tical imperfections. In Section an explicit formulus 
of the Wigner function distribution of the conditional 
State is given. In Section [VII our formahsm is appUed to 
the CW scheme, and numerical examples of the Wigner 
function and mode matching design charts are given. In 
Section rVlII the pulsed scheme is analyzed. Section IVlIII 
concludes with a few remarks. 



II. TIME- AND BAND-LIMITED SIGNALS 
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FIG. 1: Scheme of the non-Gaussian operation induced by 
the measurement by the on/off detector. 



Consider a hght beam in a single transverse mode of 
the optical field with a continuous spectrum. We denote 
the positive-frequency part of the field operator by 



à(í) = 



1 

2^ 



(1) 



Here a(wo + í^) is the annihilation operator for the Fouricr 
component at angular frequency loq + Q, where uq is 
the center angular frequency of the spectrum of a làser 
source. The operator obeys the continuum commutation 
relation 



[à{üJo + n), à^{ujQ + n')] ^ 27r5(si - n'). 



(2) 



The timc dependent field operator à{t) is defined in the 
interval (—00,00), and obeys the commutation relation 



[à{t),à^{t')]=S{t-t'). 



(3) 



Suppose that this transverse mode is excited into an 
ideal band-limited squeezed vacuum state. Let the finite 
bandwidth B in Hz, and be centered at the optical fre- 
quency ujo/2Tr. We assume that B <^ üJo/27r. Such a 
state is mathematically represented by 



Ir) =510) 



(4) 




à{üJQ + fl)à{Lü() — íl) 



-a'iuJo + Q)à''{üjQ — 51) 



(5) 



where 7(Í2) = 7(— O) is a frequency-dependent real 
squeezing parameter. This state can be obtained from 
the output of a degenerate parametric amplifier. The 
squeezing bandwidth is determined by the degree of 
phase mismatching. A cavity is often used to enhance 
the nonlinear interaction. In such a case the state Eq. 
(0)) describes an ideally simplified output state from the 
cavity provided that the parametric oscillation is oper- 
ated below threshold. The bandwidth is then given by 
the cavity bandwidth. Its precise modeling is described 
in Section IVII By the way, the band limitation becomes 



essential not for describing the squeezing process but for 
analyzing the homodyne current from detector electròn- 
ics, which has a finite, usually not so wide, bandwidth. 
It is this current that defines the measured mode. This 
point is discussed later again. 

For mathematical convenience we introduce the fol- 
lowing operator in the rotating frame about the center 
frequency ujo, 



A{t) = à{t)e"^°' = 



1 

2^ 



dnA{Q)e 



(6) 



where A{Çï) = à(a;o + íl). The squeezing operator is then 
rewritten by 



/ 1 j{n) 



-it(í7)it(_í7) 



(7) 



Using the squeezed source thus described, we now con- 
sider a scheme depicted in Fig. ^ The squeezed vac- 
uum beam is splitted via a beamsplitter with a small 
reflectance 1 — r. The reflected beam in path B is guided 
into an on/off detector. The main transmitted beam at 
path A, the signal beam, is measured by a homodyne de- 
tector. The "on" signals at the on/off detector are used 
as the triggers to select homodyne events. It is assumed 
that both the on/off and the homodyne detectors have 
the same measurement time duration T. It is this time 
interval T to define events in the time domain. Actually, 
it is not until the measured modes (states) by the on/off 
and homodyne detectors are specified that the quantum 
States of the trigger photons and the conditionally se- 
lected homodyne events are clearly defined. In other 
words, it is neither the temporal width of the input pulse 
nor its spectral bandwidth that directly determines the 
characteristics of the conditional events. 

In order to analyze this scheme we should find an ap- 
propriate orthonormal basis set for describing simulta- 
neously the squeezing, the homodyne detection process, 
and the on/off detection process. The squeezing and ho- 
modyne detection of it are described most easily in the 
frequency domain, whereas the on/off detection is based 
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on the photon counting events appearing in the time do- 
main. The continuum set of the frequency modes {^(Í7)} 
satisfying the commutation relation Eq. (O is not an ap- 
propriate set. 

Band-hmited signals to [~B/2,B/2] in Hz can be 
uniquely represented by 



ÀB{t) = ^ r dr!i(íí)e-^^*, 



1 l, 

A;— — oo 



(8) 



where is the sine function 

/— sin7rB(í - -I) 
9kit) = Vb ^ — 

satisfying the orthonormal relation 



(10) 



This set, however, is not orthogonal on a finite time dura- 
tion. So the use of this set is not convenient for describing 
the events defined on a finite time duration. 

Time-hmited signals to [-T/2, r/2], on the other 
hand, can be uniquely represented by 



Ít(í) = ÀT,k(t'k{t), 



(11) 



k— — oo 



where 



1 , 27rfcí, T T 

Mt) = < VT'''''^-'^^ (12) 



otherwise, 
satisfying the orthonormal relation 

f.T/2 



T/2 



(13) 



These functions are defined by the solutions of the inte- 
gral equation in the time domain 

, s , . í^'"^ , ,sinnB(t -t') , 

Xk{c)^k{c,t)^ dt' ^ ' ^k{c,t'), 15 

J-T/2 Tr{t - t') 

or equivalently in the frequency domain 



X/c(c)$fe(c,íí) = 



dü' 



sm 



-ífeíCÍí')- (16) 



Actually, given any T > and any i? > 0, one can find 
a countably infinite set of real functions for the integral 
equation 



(9) Xk{c)Sok{c,x) 



^ sinc(a; - y) i i / i 

dy—7-_ —Sok{c,y), \x\<l, 



-1 



(17) 



7r(x - y) 

and a set of real positive numbers satisfying 

l>Xo(c)>xi(c)>X2(c)>···, (18) 



where c = ■kBT/2. The solution Sak{c,x) is called the 
angular prolate spheroidal wave functions whose proper- 
ties are summarized in Appendixp The eigenvalues Xfe(c) 
are expressed by using a second set of solution iïp].^ (c, x), 
called the radial prolate spheroidal wave functions, as 



Xfc(c) 



2c 

TT 



(19) 



i?Q^^ (c, x) diffcr from 5oa:(c, x) only by a real scale factor. 

Using the above functions, the complete and orthonor- 
mal set on the finite bandwidth |Í7| < ttB can be given 
by 



r}k -\- 1 O 
$,.(c,í7) = \/— ^5ofe(c,— ), A: = 0,1,···. (20) 



These functions satisfy 

$fc(c,0) = (-l)^$fe(c,-r!). 



(21) 



The Fourier coefficient Àt^u is defined by 

V 1 J~T/2 



(14) 



using the operator in Eq. ©. Although this set is use- 
ful to describe the events in the time domain, it makes 
the analysis complicated when the signal is band-limited. 
Such a band-limited analysis becomes essential in de- 
scribing the homodyne current from detector electrònics 
with a finite bandwidth. 

The appropriate set to expand the time- and band- 
limited signals is known as prolate spheroidal wave func- 
tions 18]. These functions frequently apprear in diverse 
problems in physics and mathematics, and were first 
introduced to quantum òptics by Zhu and Caves |17| . 



/ dr!$fe(c,r!)$í(c,íí) = 4 



(22) 



The Fourier transforms of the functions (c, íí) give 
mode functions in the time domain. They are explicitly 



1 r" 

*fe(c,í) = TT / dmk{c,ü)e 



-int 



2í, 



= ^i2k + l)S(-z)'=iï«(c, l)Sok{c, ^)(23) 

The functions ^'^(c, í) are defined in the interval 
(— oo,cxd) and are orthonormal. 



dí*fe(c,í)4',(c,í)* = 4 



(24) 



5 



They also keep the orthogonality relation over the inter- 
val [-r/2, T/2], 



T/2 



-T/2 



dmk{c,t)-í'i{c,t)* = Xk{c)Ski. 



(25) 



For a fixed value of c the Xk (c) fall ofF to zero rapidly with 
increasing k once k has exceeded 2c/7r (= BT). A small 
value of Xfc(c) implies that í'fe(c, í) will have most of its 
weight outside the interval [—T/2, T/2] whereas a value 
of Xfe(c) near 1 implies that í'fe(c, í) will be concentrated 
largely in [—T/2, T/2]. The inverse Fourier transforma- 
tion is 



T/2 

-T/2 



iÇlt 



(26) 



The time- and band-limited fields are now quantized 
in tcrms of these mode k. We thus introduce the discrete 
set of the operators 



Ak 



27r 



-itB 



ttB 



which obeys the commutation relation 



kl- 



(27) 



(28) 



The operators À{Çl) confined in the finite bandwidth B, 
such as the ones appearing in the squeezing operator Eq. 
ÍTI), can be expanded as 



i(o) = ^ifc$fc(c,rí). 



(29) 



fc=0 



Let us assume a rectangular squeezing spectrum ^(Jl) — 
7 in Eq.(|Zj). The squeezing operator then factors into a 
product of single mode squeezing operators 



k=0 



Sk ~ exp 



with the squeezing parameter for mode k 



(30) 



(31) 



The multimode squeezed vacuum state Eq. (0J is now 
represcnted by 



to know what kind of modes are actually measured in the 
detectors. In this section we study the homodyne detec- 
tion from this point of view. The most commonly used 
homodyne detection scheme is the so-called balanced ho- 
modyne detector. Its typical scheme is shown in Fig. 
13 In this scheme the signal field A{t) is first combined 
with the LO field A^{t), via a 50:50 beam splitter and 
the two output beams are converted into photocurrents 
at each photodetector. The photocurrents from the two 
detectors are balanced to produce the difference current. 

If the photodetector has a (5-function response in time 
so that it covers the iníinite bandwidth, the photocur- 
rents are produced intantaneously. The difference cur- 
rent can then be given by 



Í{t) ^ À^\t)À{t) + À\t)À^{t). 



(33) 



In practical case, however, photodetectors themselves 
have a non-í-function response and have a finite band- 
width. Furthermore the difference current is often electri- 
cally amplified, and is finally analyzed. The bandwidth 
of homodyne detector electrònics is also limited, usually 
to a few hundred MHz at most. These effects can be mod- 
eled by the response function h(t) such that the instanta- 
neous difference current is filtered through this function 
[ï^ l2(ïl | . The operator for the filtered current is given by 



Íh{t) = / dTh{t-T)Í{T). 



(34) 



This current is integrated over the time duration 
í-T/2,T/2] 



T/2 



dtihit). 



(35) 



-T/2 



This operator directly corresponds to the final observ- 
able in the homodyne detection in the present context. 
Explícit calculation of this current Ith will be made for 
each specific physical model in later sections. 

In this section we first assume an ideal homodyne de- 
tector with the infinite bandwidth, i.e. /i(í— t) — S{t — T), 
and derive an explícit expression for the final observable. 
Since the LO field is usually a strong classical field so 
that its quantum fluctuations can be neglected. We fur- 
ther assume that the LO field has the same or narrower 
bandwidth than that of the squeezing so that the LO 
field can be band-limited at least to < ttB. One may 
then expand it as 



ir; = \rk) 
fe=o 



rk 



Sk |0) 



(32) 



III. HOMODYNE DETECTION 

In order to analyze the mode matching issue in the 
non-Gaussian operation depicted in Fig. ^ it is essential 



À^it) = ^iL^fc(c,<) = ^4e'^í'fc(c,í). 



fc=0 



(36) 



fc=0 



One should, on the other hand, note that the field op- 
erator at the signal port A{t) in Eq. can NOT be 
band-limited because the ideal homodyne detector can be 
sensitive to the modes even outside the bandwidth of the 
input squeezed beam, and these modes add the vacuum 
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FIG. 2: Scheme of balanced homodyne detection. 



fluctuations to thc homodyne currcnt. So we decompose 
the field A{t) into a sum of thc foUowing two components 



i(í) = ÀB{t) + Àv{t). 



(37) 



Here the first component ÀbÍí) is band-hmited to the 
squeezing bandwidth, and hence can be expanded as 



1 

2^ 



ttB 

TVB 



dr!i(íí)e-'"* =^ifcí'fc(c,í). (38) 



fc=0 



The second component Àv{t) consists of the remaining 
frequency modes 



Àv{t) 



1 

2^ 





rOO 


/ dn + 


/ dn 


J — oo 


JttB 



À{ü)e 



-iÜt 



(39) 



Accordingly the whole Hibert space is divided into two 
subspaces: one is Tí^ for the field modes inside the 
squeezing bandwidth |íi| < 7r_B, and the other Ti^ for the 
field modes outside the squeezing bandwidth |Í7| > nB. 
The operators Ab(í) and {Ak} act on Ti^ , while the op- 
erator Av{t) acts on . 

The ideal homodyne current is then written as 



oo oo 



oo 

+ ^4e-'^vl,*(c,í)i^(í) + h.c., 



(40) 



fe=0 



which is defined on the whole Hilbert space Tí^ ® Tí^ . 
This current is integrated over the time duration T to 
define the event signals in the time domain 

.T/2 

Ít = dti{t) 

J ~T/2 



= V2^aL^fc(c)ifee-^* 

oo 

+ V2 5]a^Vxfe(c)[l-Xfe(c)]14e-*^ 



h. c, 



(41) 



where we have used the orthogonality relation Eq. 125|l . 
and have introduced the quadrature operators on 



Vk 



1 



T/2 



\JXk{c)[l - Xfe(c)] ^-T/2 
which satisfies the commutation relation 



dtAvmi{c,t), (42) 



(43) 



Define the quadrature operators on TL^ for the modes 
inside the squeezing bandwidth |n| < -kB 

lf(0) = -i=(i,e-^^ + 4e'^), (44) 

and the ones on TC^ for the modes outside the squeezing 
bandwidth líil > ttS 



and rewrite Eq. as 

fT/2 



(45) 



It = 



T/2 
oc 



dí/(í) 



V2^4xfc(c)^f(^) 

oo 



k=0 



This integrated homodyne current operator defines the 
mode which is actually observed in the time domain anal- 
ysis, the so called LO matched mode pïL l2^. 

We define the quadrature operator for the LO matched 
mode in the foUowing way. Define the weight factors 



"fcXfc(c) 



(47) 
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and 



where 



V _ afcVxfc(c)[l-Xfc(c)] 

"^k — ~ 



(48) 



\ k=0 



such that 



k=0 



(49) 



In order to simplify the notation on tlic Hilbert space 
(g)7i^, let us further introduce vector notations 



X(<^) 



X2{4>) 



xli^) 
xfW 

x;f{^) 

xï{<l>) 



£2 



and 



x(0) 









a;o(0) 






Xi{(p) 










4(^) 
a;y(0) 



where 



and 



XU^)\xl{cb)) ^ xl{cb)\xl{cb)) 



Xn<P)\xïm^xXW\xl{<P)) 
Here we have 

[lfc(0),XK0 + |)] 



(50) 



(51) 



(52) 



(53) 
(54) 
(55) 



The integrated homodyne current Eq. (|46|l is then writ- 
ten as 



It = 



2j2Hfxkic)Xoi^) 

\ fc=0 



(56) 



^o(<^) 



X(0) 



(57) 



is the quadrature operator for the LO matched mode, 
which is defined on the entire Hilbert space Tí^ (8) Tí^ . 
It is this operator that corresponds to the final observ- 
able in the homodyme detector for the time domain non- 
Gaussian operation. In the foUowing calculation we only 
need quadrature operators at </> = 0, those are simply 
represented by 



Xk 
Pk 
Xo 
Po 



= XkiO), 

= Xkil), 

= ^o(O), 



(58) 



The POVM and the resulting measurement statistics 
for the integrated homodyne current is described by using 
the eigenstates of Xq 



Xo Wo) = Xq 14) 



(59) 



The round ket notation | ) is used to discriminate these 
eigenstates from the ones for the quadrature operators 
which describe the input squeezed {XQ.Xf, ■ ■ ■} and the 
vacuum {Xq , XY j ' ' '} fields 



Xk\xk) = Xk\xk) , (fc = 0, l,---)- 



(60) 



In order to represent the quantum state of measured 
mode, and to calculate the corresponding Wigner func- 
tion, we have to derive the formula to connect {l^g)} with 
{|xo),|a;i),···}. For this purpose we consider a real uni- 
tary transformation 



(61) 



which includes the linear relation of Eq. H57|) in the Oth 
component. We denote this equation as 







eo 


El 


£2 


Xi 






Ull 


"12 


X2 




«20 


U21 


"22 









Xi 




X2 



X = ux. 



(62) 



Such a unitary transformation is not unique, but here 
we need not to know its explícit components. As in Eq. 
(|5Ü)l. the eigenstates of Xk are denoted by 



Xk\x'k)=x'k\x'k), (fc-Cl,---)- 



(63) 



The two sets of operators {Xk} and {Xk} can be con- 
sidered to act on the same Hilbert space 
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fact, the tensor product states of all the modes for the 
two sets of eigenstates are related with each other by 



\Xk) 



because 



fc=0 

= cxp(-í*x • P) |0) 

= cxp(-í*xUt • UP) |0) 

= cxp(-í*x' • P) |0) 



fc=0 
= |X'): 



where 



x' = Ux. 



(64) 



(65) 



Given an input multimode field in a quantum state p 
on Tí^^Hy , the integrated homodyne current delivers in- 
formation only on the LO matchcd mode in p. Therefore 
what we actually observe is a reduced density operator 
on the subspace Ti}" spanned by \x'q) 



dx[ / dx'2---{x[,X2,---\p\x'i,x'2,---)- (66) 



The matrix element in terms of the quadrature eigen- 
states for the LO matched mode on the subspace is 
given by 

{x\p\y) = j dx'g J dx[ j dx'2---5{x-x'^) 
X (a;Q, x^, • • 'I p\y^x-^^,X2^ • • •) 
= /<ix'í,.-...x) 

^ i?^\p\y,x\,X2,- ■ ■) , (67) 

where we have used the relation of the eigenvalues for the 
LO matched mode 



(68) 



i=0 



We then transform the variables x' = *(xó, a;'^, Xj, • • •) 
into X = ^{xo, xi,X2t· ■ •) by Eq. Firstly we have 



dx' = 



(69) 



because the Jacobian of this variable transformation is 
|detU| = 1. Secondly 



(x'l = (x| 



(70) 



according to Eq. (|64|l . Thirdly the quadrature eigen- 
states |ï/, x'i, x'2t - ■ •) is converted as 



|ï/, Xi^X2-, 



\xk + (y - *e • x)efe) 



fc=0 

|x+ (y - 



(71) 



V 




x[ 




X2 


= u 







xo + {y- Xo)eo 
xi + {y- x'g)ei 
X2 + {y- x'ç,)e2 



(72) 



which can be obtained by eliminating the terms 
X]fc°=i '^kjx'f. in the foUowing two equations 



ut 



y 




x\ 




4 





and 



2^0 




x'q 




X\ 


= ut 


x\ 




X2 


X2 













^ly + L·k^iUkiXk 
e2y + Efeli"fc24 



eo4 + Wfco4 

£22:0 + IJfc=l Uk2Xi, 



(73) 



(74) 



The matrix element Eq. H67|l is finally given by 

(x| p\y)^ j dxS{x - *e • x) (x| p |x + (y - x)e) . (75) 

This formulus allows one to calculate the statistics of 
the LO matched mode in the time-integrated homodyne 
detection (the left-hand side), using the quantum state 
originally represented in terms of the plorate spheroidal 
wave function modes (the right-hand side). 

IV. ON/OFF DETECTOR 

In contrast to that homodyne detectors can be imple- 
mented in the near ideal condition at least for the near in- 
frared wavelengths at present using Si p-i-n photodiodes, 
on/off detectors usually suffer from imperfect efficiency 
and dark counts. In the present context where an on/off 
detector is used to select events, dark counts essentially 
influence the quality of the non-Gaussian operation. 

The number operator for photons arriving at the de- 
tector during the interval T is 

.T/2 

fiT = I dta^{t)à{t) 

J-T/2 
I.T/2 

dtÀ\t)À{t). (76) 

-T/2 

Although the signal field A{t) appearing in the homo- 
dyne current cannot be band-limited to the squeezing 
bandwidth, it CAN be here, because the vacuum field 
components outside the squeezing bandwidth does not 
induce photon counts. So A{t) can be replaced by A_b(í), 
and can be expanded as in Eq. H38|) . We then obtain 



riT 



k=0 



Xk{c)ÀlÀk. 



(77) 
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For mode k we define 

fík 



nk \nk) 



nk \nk) 



(78) 



For a photon counter with a finite quantum efficiency rjk 
for mode fc, the POVM element registeringn photocar- 
riers due to photons in mode k is given by |23l |. 



flk{n;T]k) = ^ 



77,"(l-í7fer-"|m)(m|. (79) 



The efïective quantmxi efficiency rjk for mode k here is of 
the form r]k = rjXkic), where rj is the net detection effi- 
ciency determined by the total coupling efficiency of pho- 
tons through optical components to the photondetector 
and the intrinsic quantum efficiency of the photondetec- 
tor. This is assumed to be constant over the squeezing 
bandwidth. The factor Xfe(c) represents the weight of 
mode k on the counting interval [—T/2, T/2]. Taking the 
effect of the dark counts into account, the POVM regis- 
tering n counts is then given by 



nfc(n;ryfe,í/fc) 



(n-n')\ 



-nk{n':,r,k), (80) 



where Vk is the mean number of the dark counts for pho- 
tons in mode k. (The dark counts may occur regardless 
of the bandwidth of the input signal field, however, they 
can be taken into account by adjusting the vàlues Vk for 
the band-limited modes such that the total dark count 
rate in practical detecotrs is properly modeled.) These 
elements satisfy the completeness relation for each mode 
k, 



'^Ilk{n;rjk,iyk) = ík- 

ji=0 



(81) 



The multimode on/off detector placed at path B in Fig. 
^is finally modeled by a binary POVM with parameters 
V = {m,Vi, ■ ■ •) and = {i^o,!^!, ■ ■ ■), 



V. 



nBj(r,,i.) = (g)ni^(0;77fc,^fe), 

fe=0 



WIGNER FUNCTION OF THE 
NON-GAUSSIAN STATE 



(82) 



We now explicitly calculate the non-Gaussian output 
State conditionally obtained by the on/off detector spec- 
ified in the previous section, and derive an expression of 
the Wigner function measured by the ideal homodyne de- 
tector described in Section [llll For convenience of math- 
ematical handling, we express the squeezed vacuum state 



rk 



rk 



Sk |0) , 



(83) 



in terms of the coherent state basis \ak) on the Hibert 
space Ti,^ . The component of mode k can be expanded 
as 



{l-Xk')Í 



(fak exp 



, (84) 



with Xk — tanhrfe. As in Fig. ^ the squeezed vacuum 
state is beamsplitted into path A and B, resulting in a 

state 



{i-Xk')i r ,2 

' d ak exp 



fc=0 



where r is the transmittance. 

The photon subtracted non-Gaussian state at path A 
is given by 



pi 



1 



Tr 



Pdct 

where 

-Pdct = Tr AB 



B 



^')(AB)(^V^®n:fjr,,«.)l , (86) 



/'')(AB)(^'l^^®nBjr,,i.)' 



(87) 



is the probability of having the "on" signals. Using Eqs. 
(|S5|l and (|S3 . p\ can be represented as 



pi 



1 



with 

RkiVk, Vk) 



liïfc(0,0) 



\RkiVk,Vk) 



k=0 



(fak / (ff3kexp 



2-T 



X i^\ak\' + IPkl' j + {1 - t){1 - Tjk)P*kak 
Vk \V^ak)j^{V^í3k\. (89) 



The state is finally measured by the time-integrated 
homodyne detection to construct the Wigner function. 
As explained in Section IIIII the ideal homodyne detec- 
tor is sensitive not only to the state but also to the 
vacuum states |0^^ outside the squeezing bandwidth. So 
the input state into the ideal homodyne detector must be 



/5a = pÍ<^|0^)(0VI 



(90) 



k=0 



The time-integrated homodyne current delivers Informa- 
tion only on the LO matched mode, whose state is rep- 
resented by the reduced density operator pA defined by 
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Eq. H66|) . Its matrix element in terms of the quadrature 
eigenstates is given by Eq. l(7S|) . The Wigner function of 
the reduced density operator for the LO matched mode 
is then given by 



Wix,p) 



1 

2^ 



PA 



(91) 



where the matrix element is given by 











(-1 


PA 


X — 


1) 



dxS{x + - - *e • x) (x| PA |x - ^e) 

= j j d^^5{x + I - . x^ - • x^) 
X (x^l Pa - ^w^) 

x(xV|0V)(0V|xV-^wV). (92) 

In the second equality we have replaced the abbreviated 
notations {e^} and {xk} introduced in Eqs. (|ïïT|l and ifS^ 
with the explícit ones {w^,w^} and {x\,x^}. Substitut- 
ing the expression of Eq. (|88|l into the above equation, 
we have 











(-1 


PA 


X — 


i) 



oo „oo 

i^„T^- / . V 



oo 



n 

/=0 



dxiexp{-iqwYxY) 



x{xr\or){Qr\xr-^wr) 

Y[ / dxlexpi-iqwfxl) 
x{xl\RkiO,0)\xl-Íwl) 
Yl / dxlexpi'-iqwlxl) 



k=0 



(93) 



where we have used the Fourier expansion for the delta 
function in Eq. H92|l . After straightforward calculations 
of Gaussian integrations, the Wigner function can be fi- 
nally represented as 



W{x,p) 



1 



Pde 



[R{x,p;0,0)-R{x,p;v,'^)], (94) 



where 



R{x,p;t],u) = 



X exp 



C-iv) C+iv) 



(95) 



with 



C±(r?) = l±r 



2Ai 



7±(í?fc) = 1 T Afe ± (1 - T)r]kXk, 



and 



k=0 



1+ivkh-ivk) 



(96) 



(97) 



(98) 



and 



R{x,p;0,0) 




air 1 



(1 - T + a'ÍT){l - T + (t\t) 



(99) 



with 



oo 

k=0 



(100) 



The probability of having the "on" signal is given by 

Pdct-l-AA(r7,iy). (101) 

Thus the relevant quantum states can be represented 
in terms of the discrete set of the prolate spheriodal 
wave function modes. This discretization becomes more 
prominent as the time- and band-limitation gets more 
stringent. In fact, for BT ^ 1, only the Oth mode be- 
comes dominant as shown later. If one could make the 
photon-subtracted squeezed state pure in this Oth single 
mode, and could selectively measure the Oth mode by the 
homodyne detector, then the ideal mode matching should 
be realized. But for this, it is not sufíicient to simply put 
BT <C 1. Some additional care must be taken, depend- 
ing on specific physical models actually used. So in the 
foUowing two sections, we apply the result in this section 
to two kinds of models: CW and pulsed schemes. 

VI. CW SCHEME 

A typical example of CW scheme is an optical paramet- 
ric oscillator cavity containing a x*-^-* nonlinear medium, 
continuously pumped by a single mode field at optical 
frequency ^uj^. A pump photon is down-converted into 
quantum correlated twin photons at optical freqüències 
LúQ -\- Çl and — í^, resulting in an frequency entangled 
state of squeezed vacuum via cavity enhancement. In 
a commonly used temperature-controUed phase match- 
ing scheme, such twin photons are generated over a wide 
range of freqüències something like a few tens of GHz 
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«0 - fflpSR 



(O 



ft>o+;rB ft> 



FIG. 3: Frequency spectrum of the squeezed state generated by an optical parametric oscillator cavity. 



around the degenerate frequency loq . A typical frequency 
spectrum is shown in Fig. |21 It consists of multi resonant 
peaks separated by the frequency o^fsr of a free spec- 
trum range. For a commmon bow-tie configuration with 
a round-trip length of ~500mm, üjpgri is about 600 MHz, 
and a width of each resonant peak is B ~ 10 MHz. Such 
a wide band CW squeezed state is a source state for the 
non-Gaussian operation here. The events in time domain 
are defined by imposing the finite time duration T on the 
on/ofï detector. 

Let us íirst consider time-integrated homodyne detec- 
tion with the infinite detection bandwidth as shown in 
Fig. ^ The LO is a CW single mode field at optical 
frequency lüq. In the rotating frame, it is written as 



a e 



(102) 



A CW signal beam is combined with the LO field, pro- 
ducing the CW current. 



1 f°° 
I{t) = — / dül{n)e 
27r J„oo 



iQ.t 



with 



I{ü) = 



27r 



À{n)e 



— Í(j) 



,icti 



(103) 



(104) 



and this current is integrated over [—T/2, r/2] 

L coo 



277 



dü 



Ain)e 



sm 



(105) 

As seen from this expression, it depends on the choice of 
T what frequency modes dominate in It- For T ^ í^fsR' 



CW single mode LO field 



BS 



I 



Signal beam '^-'«t ^^^^^ photodiode 
' T=0.5 ' 



i{t) 




FIG. 4: Time-integrated homodyne detection in a CW set- 
ting. The detection bandwidth is assumed to be infinite. 



contribution from the modes in the resonant peaks at 
Lüo ± nwpsR {n = 1, 2, ...) is negligible. For T > B^^, the 
modes inside the resonant peak centered at luq dominate. 
For an intermediate region, Wpgp, < T < , the modes 
in the free spectrum range, which are in vacuum states, 
also contribute to It in addition to the ones in the center 
resonant peak. Thus under the assumption T ^ 'j-'fsr' 
it is enough, concerning the time-integrated homodyne 
detection, to consider a squeezed state whose spectrum 
is confined to the center resonant peak shaded in Fig. O 
One should, however, note that this is not true for 
the on/off detector. The on/off detector is sensitive to 
all the modes regardless of the measurement duration 
T. Therefore for a good frequency mode matching, the 
tapped-off' beam must be optically filtered in front of the 
on/off detector such that only the center resonant peak 
shaded in Fig. |2|is guided into the on/off detector. This 
can be made by employing an appropriate set of filtering 
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cavities. This is actually the assumption under the model 
of band-limited squeezing introduced in Eq. © or Q. 
We further simplify the squeezing spectrum so as to be 
the flat band over [—ttB, ttB], as shown in the lower part 
of Fig. 121 which is actually the model of Eq. l|ÏÏÜ|l . 



LO matched mode measured by homodyne 
detector with the infinite bandwidth 



Before analyzing the non-Gaussian operation, let us 
consider a simple case where the squeezed vacuum state is 
directly measured by the homodyne detector without the 
photon-subtraction by the on/off detector. The Wigner 
function in this case is given simply by R{x,p; 0,0) in 
Eq. ^ and ifïÜÜjl with r = 1. Using Eq. the kth 

coefScient of the CW single mode LO field in Eq. H1U2|) 
is given by 



al 



— $fc c,0 = — 



2k + l 



B 



where Pk{x) is the kVa Legendre polynomial and 
. (fc-1)!! 



Pfc(O) 




fcü 



(fc : even), 
(fc : odd). 



The variance is then expressed as 



where 



E 



fc:cvcn 



^ {2k + l)Pu{0fxk{c) 



(106) 



(107) 



(108) 



(109) 



and 



E 

A::cvcn 

^ {2k+i)pm\um-xk{c)] 



fc:cvcn 



^ {2k+l)Pk{Qfxk{c) 

fc:cvcn 



(110) 



In the case of BT (= 2c/ -k) > 1, the first TV - BT 
modes have eigenvalues Xk{c) near 1. These modes arc 
foUowed by a transition region of /C ~ 2/7r^ ln(27ri?T) 
modes, in which the eigenvalues fall from near 1 to near 
0, as shown in Fig. [S] Beyond the transition region, the 
remaining modes have eigenvalues that are very close to 
0. This eigenvalue spectrum and the above Eqs. H109|) 
and (|110|l mean that w^ ~ 1 and w^ ~ 0. Thus in 
this limit (a long enough measurement interval T) the 
original squeezing characteristics can be directly observed 
without any degradation due to the vacuum modes. 



A^~4ln(27ieT0 




FIG. 5: Eigenvalu spectrum of Xfc(c) for BT 3> 1. 



k 


BT = 0.1 


BT = 0.5 


BT = 1.0 


BT = 3.0 





0.09973 


0.46780 


0.78340 


0.99890 


1 


0.00027 


0.03183 


0.20502 


0.96869 


2 


0.00000 


0.00037 


0.01136 


0.73284 


3 




0.00000 


0.00021 


0.26248 


4 






0.00000 


0.03478 


5 








0.00221 


6 








0.00009 



TABLE I: The first few eigenvalues Xk{c) for severa! vàlues 
of BT. 



In the case of BT ^ 1, on the other hand, the Oth 
eigenvalue scales as Xo(c) ^ BT and Xo(c) ^ Xfe(c) for 
A: > 1, namely, only the Oth mode contribuïes with the 
weight BT . The first few eigenvalues Xfc(c) for several 
vàlues of BT are given in Table ^ which is borrowed 
from Zhu and Caves (Tt'I. This means that in Eq. (|1Ü8|I . 
w^ ~ BT < 1 and w^ ~ l-BT w 1. The second term in 
Eq. p()8|l represents the vacuum fluctuations due to the 
modes outside the squeezing bandwidth. Actually if one 
measures the squeezed state in an interval much shorter 
than ~ B~^ , one will observe wide range of frequency 
modes in vacuum states. When this term becomes dom- 
inant, the squeezing characteristics cannot be observed 
clearly, being covered by the vacuum fluctuation noise. 
This could be a serious restriction on optimization of the 
mode matching. In the next subsection we consider a 
scheme to remove this restriction by using an electrical 
low pass filter. 



B. LO matched mode measured by homodyne 
detector with a flnite bandwidth 



Let us consider a scheme shown in Fig. [ü] where an 
electrical low pass filter is inserted into the scheme of Fig. 
2] The intantaneous homodyne current Eq. H33() is first 
filtered by the low pass filter with the same bandwidth 
B matched to the squeezing spectrum, and is then inte- 
grated over the interval [—T/2, r/2]. The filtered current 
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CW single mode LO field 
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FIG. 6: Time-integrated homodyne detection with a low pass filter in a CW setting. 
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FIG. 7: CW scheme of tl·ie non-Gaussian operation by pl·ioton subtraction with the on/ofï detector. 



can be simply represented as 



lBÍt) = 



47r2 



iílt 



It is finally time-integrated over the interval T 



(111) 



Itb = 



27r 



ttB 



(112) 

In this equation the field A{Çl) is band-limited and hence 
can be expanded in terms of Ak by Eq. I|29() , resulting in 



Itb — 



k=() 



ttB 



Í2T 



ttB 



(113) 



where the definition Eq. ll^ is used by setting = 0. 
By using the relation Eq. 116|l . the above equation is 
rewritten as 

L 

íTB = ^J2^kXkmk{c,o). (114) 



So in this case the observed LO matched mode is specified 
only by the quadrature operator on the Hilbert space Ti.^ 



(115) 



k=0 



where 



e/c 



Xfe(c)$fc(c,0) 



J2xk{c)Hk{c,oy 

\ k=0 



(116) 



which is non-zero only for even fc, because of Eq. H106() 
and Eq. pU7|l . The observed variance Eq. H1U8|) thcn 
reduces to 



2 -2t 



(117) 



which is independent of the intcgration timc T. Thus 
in this case one can observe the intrinsic characteristics 
of the squeeezing regardless of the choice of T. This is 
essential for achieving stringent mode matching in the 
regime of BT < 1. 



C. Numerical examples and mode matching design 

chart 

After all, we consider the CW scheme depicted in Fig. 
[3 The CW squeezed beam is splitted with reflectance 
1 — r = 0.1. The 10% of the squeezed beam is tapped off, 
and is guided into the on/off detector, which is opened for 
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the time duration [—T/2, r/2] by electrical gating. This 
defines the discrete events in the time domain. The aver- 
age interval between the successive trigger events ("on" 
counts) is assumed to be long enough compared with T. 
In the homodyne detector, the CW signal bcam is com- 
bined with the CW single mode LO field, producing the 
CW current /(í), and this is filtered by the low pass fil- 
ter matched to the squeezing bandwidth. Only when the 
"on" signal is sent from the trigger channel, the filtered 
current /^(í) is integrated over [—T/2, T/2] synchronized 
with the gating signal. The quadrature operator of Eq. 
(IllSfl based on the integrated current Itb is used to con- 
struct the Wigner function of the conditional statistics. 

As a typical model of squeezing, we assume 7 = 0.35 
(3 dB squeezing) and B = 10 MHz. The net detection 
efSciency of the on/off detector 77, which appears in the 
effective quantum efïiciency for mode k as rjk — ?7Xfc(c), 
is determined by the total coupling efficiency of pho- 
tons through filters and couplers to the photondetector, 
and the intrinsic quantum efliciency of the photondetec- 
tor. So this can be a small value something like <0.5. 
The mean number of dark counts Vk for mode k can 
be converted into the dark count rate n[counts/s] by 
J2k'^k = nT. We vary this n [counts/s] to evaluate the 
dark count effect. 

In Figs. ISltroughim we show the vàlues of the Wigner 
function at the phase space origin as a function of the 
dark count rate n [counts/s] for several vàlues of BT. 
The mode weights {xfc(c)} are given in TableQ] BT = 
means the single mode case with perfect mode matching. 
Figs. IHl El El and ^] corresponds to the different de- 
tection efficicncies of the on/off detector, 77 = 0.01, 0.1, 
0.7 and 1. From these figures, we can know the thresh- 
old for the on/off detector dark counts below which one 
can expect to observe the negative dip of the Wigner 
function, which is a sign of the non-classicality of the 
non-Gaussian output state. As the detection efïiciency 77 
becomes smaller, the threshold dark count for the nega- 
tive dip also gets smaller. 

For larger vàlues of BT, the mismatch between the 
LO matched mode and the photon mode observed by the 
on/off detector becomes more serious. For BT — 3.0, one 
could not attain the negative dip of the Wigner function 
even with the ideal on/ofï detector. Actually, in the case 
of BT ^ 1, the first Af ^ BT modes have eigenvalues 
Xk (c) near 1 , and then the eigenvalues fall from near 1 to 
near rapidly as k increases, as shown in Fig. |S1 All the 
first Af ^ BT modes cause trigger signals at the on/ofï 
detector, which makes the conditional state at the signal 
port a highly mixed state, because the on/ofï detector 
cannot discriminate which mode a photon comes ïrom. 
On the other hand, the homodyne detector sees only the 
LO matched mode, which is a particular combination 
of the first Af BT modes, and provides the homodyne 
statistics ïor the single mode quadrature operator defined 
byEq. (tTTÏÏll . 

In the case of BT <C 1, on the other hand, Xo(c) ~ BT 
and Xo(c) ^ Xfc(c) for fc > 1. So only the Oth mode is 
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FIG. 8: Vàlues of the Wigner functions at the origin of the 
phase-space, versus the dark count rate n [counts/s], under 
the condition of 3 dB squeezing, 77 = 0.01, and r = 0.9 
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FIG. 9: Vàlues of the Wigner functions at the origin of the 
phase-space, versus the dark count rate 7i [counts/s], under 
the condition of 3 dB squeezing, 77 — 0.1, and r = 0.9 



dominant both in the trigger channel and the homodyne 
detector: more precisely, the POVM element for the "on" 
signal U^^{t],L·>) in Eq. and the quadrature eigen- 

state \x) on the Hibert space describing the homo- 
dyne detector (see the text from Eq. l|ÏÏÜ|l to Eq. ifTsl) '). If 
one further takes a small reflectance for the tapping beam 
splitter, making the probability of detecting more than 
two photons at the on/off detector very small, then the 
trigger photons are projected onto an almost pure sin- 
gle photon state. The homodyne detector measures the 
same pure state with almost perfect efhciency, attaining 
the best mode matching. The small tapping fraction, 
combined with a small value of Xo(c) ^ BT, means a 
small efficiency for the on/off detector. This is usually 
unwanted, but in the present context this simply results 
in the reduction of the number of selected events. If the 
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FIG. 10: Vàlues of the Wigner functions at the origin of the 
phase-space, versus the dark count rate n[counts/s], under 
the condition of 3 dB squeezing, rj = 0.7, and r = 0.9 
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FIG. 11: Vàlues of the Wigner functions at the origin of the 
phase-space, versus the dark count rate n[counts/s], under 
the condition of 3 dB squeezing, rj = 1, and r = 0.9 



number of true trigger events can be kept larger than 
the dark counts, this reduction would be an acceptable 
sacriface to attain better mode matching. 

One should note that it is essential to filter the ho- 
modyne current by the low pass filter matched to the 
squeezing bandwidth B before integrating over the in- 
terval [—T/2, r/2]. If the ideal homodyne detector with 
wider enough bandwidth than B would be used, the vac- 
uum fluctuation would also become dominant in the ho- 
modyne statistics especially for shorter T. This compo- 
nent does not have any mode overlap with the trigger 
photons at the on/off detector, and would degrade the 
mode matching quality. 

Figs. El El and 1 1 41 represent the Wigner function dis- 
tributions at the point of dark count n — 500 [counts/s] 
in Fig. M{V = 0.1) for the three kinds of BT = 0.0, 0.5, 




FIG. 12: Wigner function distribution for BT = 0, i.e. 
the ideal single mode case, at the point of dark count n = 
500 [counts/s] in Fig. |Ü1(3 dB squeezing, r = 0.9, r; = 0.1) 




FIG. 13: Wigner function distribution for BT = 0.5 at the 
point of dark count n — 500 [counts/s] in Fig. |Ü1(3 dB squeez- 
ing, r = 0.9, 77 = 0.1) 



and 1.0, where one can still expect the negative dip. 



VII. PULSED SCHEME 

In this section we consider a scheme using Fourier- 
transform limited short pulses. A pulse width of the 
fundamental field (LO field) is typically 1 ps or shorter, 
so that the spectral width is about 1 THz or wider. 
The pump field is the frequency-doubled LO with al- 
most the same spectral width. A nonlinear crystal is 
pumped in a traveling-wave configuration in a single 
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FIG. 14: Wigner function distribution for BT = 1.0 at the 
point of dark count n — 500 [counts/s] in Fig. |Ü|(3 dB squeez- 
ing, r = 0.9, Ti = 0.1) 



path. Then unlike the CW configuration, the down- 
converted (squeezed) field has generally difFerent mode 
properties from the ones of the pump and the LO fields. 
The finiteness of the nonhnear crystal broadens the 
phase-matching condition, producing a wider bandwidth 
of squeezing than that of the LO. The group velocity 
mismatch increases the pulse duration, and makes the 
field apart from the Fourier-transform hmit. The group 
velocity dispersion induces a frequency charp. These dis- 
crepàncies between the signal and the LO field proper- 
ties degrade the degree of mode matching. In order to 
overcome this problem, the usage of optical filters with 
narrow bandwidth in the trigger channel was enphasized 
in 0, 0, 0, Q|. In these works, a biphoton source 
State 



Ip)ab = |0)aIO)b 



cfkdujcfk'duj' 



(118) 



is assmTied. Here |k, íí;)^ means the singie photon statc 
with momentum k and frequency to in path j (=A, B), 
and the function í'(k, uj, k', uj') characterizes the correla- 
tion properties of the biphoton state. Since we are con- 
cerned here with frequency mode matching, we shall sup- 
press the momentum representation. The trigger photon 
at path B is sclccted by frequency fihers, and is then 
detcctcd by an ideal photon counter. When the pho- 
ton countcr clicks, the trigger state is projected onto a 
POVM element 



Òon = / dujF{uj) . 



(119) 



where F{uj) includes filter characteristics and the effec- 
tive quantum efhciency of the photon counter for each 



frequency mode. The conditional signal state at path A 
is then given by 



PA 



TrB (tíon |p)(AB) {P\ 
Tl-AB (rïon \P)(AB) {p\ 



(120) 



In lïjl, the correlation function 

r{u,iü')^Tr[pAà\Lü)àiu)], (121) 
and the modal purity 

P = Tr(pi), (122) 

were introduced, and a measure of mode matching M 
is defined using the correlation function. The quality of 
mode matching was analyzed in terms of these M and 
P. However, it has been unclear how to obtain explícit 
expressions of the quantum states when these quantities 
are given. 

Our interest here is a more direct quantity than M and 
P, that is, the Wigner function of pA finally obtained by 
the homodync detection. In addition, the source state 
is not the photon pair but the wideband squeezed state, 
which is beyond the first order perturbation theory. The 
source state is a flat band squeezed state |ro) = Sq |0) 



dfl 



Bs 



A{Q)A{-Q) 



(123) 



where the squeezing bandwidth Bg is usually much wider 
than that of the LO B, i.e. Bg ^ B. Here we have omit- 
ted the group velocity dispersion and mismatch, because 
this simply causes slight quantitative modification for the 
filtered state, and is not essential in considering the fre- 
quency mode matching here. 

We should first derive the quadrature eigenstate \x) to 
describe the POVM of the homodync detector. For sim- 
plicity, we assume the LO field with a flat band spectrum 
over the bandwidth B 







otherwise. 



(124) 



The instantaneous homodync current is then given by 



Í{t) = ^ / rfí7/(íí)e-'"* 



(125) 



ith 



ZTT 



dü' 



+AH-ü + n')A^{ü') 



L r-^B 



ttB 



À{ü + íí')e"'"^ 



(126) 
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Note that the field operator at the signal port A(íí) itself 
is not band-limited. If the homodyne detector would have 
the infinite bandwidth, then the instantaneous current 
I{t) provides the final observable. In practice, however, 
this is not the case. Rather, the efi'ective bandwidth of 
the homodyne detector including the photodiodes and 
the electrònics, Bn is much narrower than that of the 
LO. Typically, í?h is a few hundred MHz at most while 
B >lTHz. So one actually measures the low frequency 
components within the homodyne detector bandwidth. 
We shall approximately model it as 



— / dni{fl)e 



a 

4^ 



düe-'^* / dn' 

-ttBh J-ttB 



X 



sin TiB^ít 

2^ 



X / dü 



(127) 



The field operator À{íl) is now band-limited to B, and 
can be expanded by the discrete set of the prolate 
spheroidal wave function basis. Using the formula in Sec- 
tionm we have 

/BH(í)-V2aL^^^f:x|(^)*,(c,0), (128) 

k=0 

where the quadrature operator X^{(j>) is defined by Eq. 
(|44|l . The time dependent factor means that a short 
optical pulse with a duration T ^ B~^, say ~ 1 ps, 
is converted into an electrical pulse with a duration 
Th ^ -Bfj^ ~ 0.1 ns. As seen from this equation it is 
sufficient to sample the peak value of the electrical pulse 



/bh(0) « V2aLi3H^i-S(^),í,^(c,o), 



(129) 



fc=0 



for the quadrature vàlues. The final observable is then 
given by 



^0 - y^efc^fc. 



fc=0 



where 



í'fe(c,0) 



^|*fe(c,0)| = 

\ fe=0 



(-l)V(2^· + l)Xfc(c)Pfc(c,0) 



\ k=0 



which is non-zero only for even k. The POVM of the 
pulsed homodyne detector is represented by the projec- 
tion \x) {x\ onto the eigenstates of Xq, tracing out all the 
other modes. 

If the trigger beam is projected onto the single pho- 
ton State |1) on the subspace spanned by \x), the 
ideal photon-subtracted non-Gaussian state would be ob- 
tained. Toward this limit, the trigger beam must be 
band-limited to B, since the homodyne observable con- 
sists of the frequency modes within the LO bandwidth 
|íl| < 7r_B. It can be made by interference band pass 
filters typically with 0.1-1 nm spectral width. The band- 
limited trigger beam is finally detected by the on/off de- 
tector. The on/off detector responds much more slowly, 
say in a scale of Td ^ 300 ps, than the pulse width 
T < 1 ps. But if we assume that the average interval be- 
tween "on" signals is still much longer than this Td , and 
also that the photoelectric conversion does not signifi- 
catly destroy the Fourier-transform limited pulses. Then 
the on/off detector can be represented by the state basis 
of the prolate spheroidal wave functions for the BT, as 
described in Section IÏVl 

The photon subtracted non-Gaussian state at path A 
is given by 



Tr 



PA = 



|P)(AB) {p\Í^<S>fl^M'^) 



Tr 



AB 



|P)(AB) (p|/A®n?n(^,I^) 



(132) 



Since IIqjj detects only the modes in |Í7| < ttB, and these 
modes are not quantum correlated with the ones outside 
of it, the state |p)ab '^^^ be regarded as the beam splitted 
state from the squeezed state |r) — S \ 0) where 



S 



exp I — / dQ 



A{n)A{~n) 



I exp 



k=0 



{Àl 



-A 



12 



(133) 



instead of the squeezed state |ro) = S'ò |0) of Eq. (|123|l . 
Thus the photon subtracted state pA is represented in 

(130) terms of {Ak} or {^f} and its eigenstates {|a^^)}. Its 
Wigner function can then be calculated by the formulus 
(|94|l . The mode matching consideration proceeds just 
as the CW case in Section FVI Cl Good mode matching 
is achieved when the weight coefficient eg dominates the 
other ones ei,e2,.... This can be made by setting the 
valuse BT as small as possible. In the pulsed scheme, 
however, the T is automatically set by the làser pulse 
width. It is actually impossible at present to realize 

(131) ^ shorter measurement interval because electrical gat- 
ing cannot reach such a time scale. Therefore the value 
of BT is lower bounded by the Fourier-transform limit, 
which is about BT w 1 or slightly less. Then a few lower 
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modes can contribute, as shown in Table |IJ This is al- 
ready enough to observe the negative dip in the Wigner 
function, provided that mode matching with respect to 
the other degrees of freedom can be made perfect. 

In experimental practice of the pulsed scheme, how- 
ever, it is generaUy not easy to reahze a high quahty 
spatiotemporal mode matching unhke the CW scheme 
based on cavity OPO systems. For meaningful numer- 
ical evaluations, consideration on imperfections due to 
the spatial mode mismatch should be involved. But this 
is beyond the scope of this paper. 



VIII. CONCLUSION 

We have developed a theory to design the frequency 
mode matching in the photon-subtracting operation 
by the on/off- type photon detector on the wideband 
squeezed state. It is essential to represent the POVMs 
of the on/off photon detection and the homodyne detec- 
tion in terms of an appropriate basis under time- and 
band-hmitation on the fields. Such a set is based on the 
prolate spheroidal wave functions which are complete and 
orthonormal under the time- and band-hmitation. The 
POVMs thus represented define the measured modes, 
and hence the quantum states of the trigger photons and 
the conditionaUy selected homodyne events. The mode 
matching is pursued for those quantum states. 

Our theory has been apphed to the CW and pulsed 
schemes. In both schemes, the desired condition for good 
frequency mode mathching is BT ^1. In the prcvious 
works [ÏO, IH, [Í3) Hjl ) only narrowband spectral filtering 
is emphasized. But the quantity to be set smaller is the 
product BT. In addition, the modes of interest are not 
the ones specified by the plane wave basis [10111^113.11^ 
but rather the ones specified by the prolate spheroidal 
wave functions under the time- and band-limitation. In 
the regime BT ^ 1, the discreteness of modes becomes 
prominant, and only a few lower modes of the prolate 
spheroidal wave functions are excited. For BT < 0.5, 
only the Oth mode becomes dominant with the mode 
weight Xo(c) ~ BT as seen in Table U Then the trig- 
ger channel selects a pure single photon state in mode 
if one takes a small reflectance 1 — t of the tapping beam 
splitter, and could suppress dark counts. The homodyne 
detector measures the same pure state with almost per- 
fect efficiency, attaining the ideal mode matching. 

The increase of BT and the dark counts makes the con- 
ditional state more mixed over many modes of fc, while 
the homodyne detector still projects the conditional state 
onto a pure eigenstate of the LO matched quadrature 
operator, which is a particular combination of the first 
M ~ BT modes. Then the non-classical effect in homo- 
dyne statistics will be smeared out. 

In the CW scheme where all the beams are CW, such 
as based on cavity OPO systems, the effective bandwidth 
can be set in a range of i? w 10 MHz by using an appro- 
priate set of filtering cavities placed in fromt of the on/off 



detector, and an electrical low pass filter in the homodyne 
detector spectrally matched to the squeezing. Then the 
time duration T for the desired condition BT < 1 is of 
order of sub /xs, which can readily be achieved by current 
electrical gating. A smaller T allows one more stringent 
frequency mode matching. The lower bound for T may 
be set by the temporal resolution of the homodyne de- 
tector, which is currently of order of a few ns. So the 
frequency mode matching in a regime BT ^0.1 can be 
attained in principle. The spatial mode matching can 
also be fulfiUed by carefuUy locking the cavities. The 
temporal mode matching in the sub /is is not a problem. 
Then the numerical results presented in Section IVI Cl 
can be practical design charts. One may expect to real- 
ize the negative dip of the Wigner function at the phase 
space origin for practical experiental parameters. All in 
all, the CW scheme will provide a good test-bed for the 
non-Gaussian operations based on photon counting and 
homodyning. 

In the pulsed scheme, the time duration T is auto- 
matically set by the làser pulse width, typically ps or- 
der. The band limitation comes from the LO bandwidth, 
sa ITHz, corresponding to the Fourier-transform limit 
0.5 < BT ^ 1. The time response of the homodyne de- 
tector is much slower than the LO pulse width. Then by 
sampling the peak value of the homodyne current pulse, 
one finally observes the quadrature composed linearly of 
a few lower modes of k. This makes the situation very 
similar to the CW case, and the numerical examples in 
Section IVI Cl can apply. So the frequency mode match- 
ing needed for observing the negative dip of the Wigner 
function will be possible. In practical experiments, such 
as ref. 0], however, the observed dip of the Wigner func- 
tion is still positive. This may be partly attributed to 
the imperfect spatial mode mathing between the trigger 
photons and the LO field. It is generally more difficult, 
compared with the cavity CW scheme, to control the 
wave front of the pulsed field, and to spatially match the 
relevant modes. 

Now we would like to mention the relation between our 
result and the one obtained by Grosshans and Grangier 
[ï^. They predicted the lower degree of mode matching 
for the CW scheme than for the pulsed one. This seems 
to be contradict with our result. In however, some 
points are missing. Firstly, the POVM element in the 
trigger channel was modeled by a projection onto a pure 
state 

dn\n) / dn'{íi'\, (134) 

-ttB J-ttB 

instead of Eq. H119() . But photon counters cannnot be 
sensitive to the relative phases of different frequency com- 
ponents. So the modeling by the above equation is un- 
likely, as also pointed out by Aichele et al. [Í3 . Secondly, 
the parameter X = Su>íST/2tt in |l2l |. which corresponds 
to BT in this paper, could not be arbitrary small for 
the Fourier-transform limited pulses. According to these 
points, the degree of the mode overlap ricf[{X) should be 
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reconsidered in the pulsed regime. Finally, the analysis 
on the "chopped" CW scheme did not take the time- and 
band-hmitation into acount. In the single parameter 



VcsiBT) 



1 

BT 



dü- 



ttVL 



(135) 



was used for quantifying the degree of mode matching. 
But this cannot be sufficient for treating the time- and 
band-hmited signals. More precisely the mode matching 
must be characterized by the set of mode weights 



dm^ic, n)- 



ttB 



X/c(c)$fe(c,0). 



(136) 



Therefore the upper bound 7?cff(^) < 0.825 for the CW 
scheme derived in 12] will not be the true hmit. Rather 
the more stringent matching wiU be possible as shown in 
this paper. 

Concerning to the CW scheme, we should also mention 
the effect of chopping the beams for the dm'ation T. The 
scheme in this paper is to use the CW single-mode LO 
and then to integrate the CW homodyne current over 
T. This is equivalent to using the chopped LO with the 
duration T, if the detector bandwidth is infinite. Re- 
member that the homodyne detector with the infinite 
bandwidth suffers from the vacuum fluctuations outside 
the squeezing bandwidth for shorter T. So we have con- 
sidered installing the matched low pass filter. In such a 
band-limited case, using the chopped LO will not gen- 
erally be equivalent to integrating the CW current over 
T afterward. In fact, chopping the CW single-mode LO 
modifies its spectrum from S{íl) to sm{ílT/2)/Trn. Then 
the homodyne current is also modulated accordingly, and 
is then filtered. According to the present analysis, the use 
of the chopped LO does not seem to bring any particu- 
lar mèrit. Therefore it will be enough to simply use the 
CW single-mode LO, and to electrically gate the on/of 
detetor to define the events in the time domain. 
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APPENDIX: PROLATE SPHEROIDAL WAVE 
FUNCTIONS 

In this appendix bàsic properties of the prolate 
spheroidal wave functions are sumarized for reader's con- 
venience. The differential equation 



(1 



dx-^ 



2x^ + l·j.-c'^x'^)u = 0, (A.l) 
dx 



has continuous solutions in the closed set x interval 
[—1. 1] only for certain discrete real positive vàlues, 



< /io(c) < Mi(c) < M2(c) < 



(A.2) 



Corresponding to each eigenvalue fik (c) , there is a unique 
solution Sok (c, x) such that it reduces to the A:th Leg- 
endre polynomial Pk{x) uniformly in [—1,1] as c 0. 
The functions Sok (c, x) are called the angular prolate 
spheroidal wave functions. They are real for real x, con- 
tinuous functions of c for c > 0, and orthogonal in (—1, 1). 
S'ofe(c, x) has exactly k zeros in (—1, 1), and even and odd 
according as k is even and odd. 

Alternatively 5'ofe(c, x) is also a solution of the integral 
equation 



Xkic)Sok{c,x) 



smc{x - y) 11/1 

dy — —Sok{c,y), \x\<l. 

-1 T^{x-y) 

(A.3) 

The eigenvalues Xk (c) are expressed by using a second set 
of solution R^i}{c,x) for (|A.ip . called the radial prolate 
spheroidal wave functions, as 



2c 



(A.4) 



R^^k í*^' ^) diffcr from S^kic, x) only by a real scale factor. 
From Eq. lfO|) . we have 

Xfc(c)*fe(c,í) = 

and 

Xfe(c)4'/c(c, n) 



dt' ^ ^fc(c,tO, (A.5) 

-T/2 '^[t - t') 



TvB ■ (n-n')T 

2 -Mc,í^')- (A.6) 



dn'- , 



The orthogonal condition of Eq. H25(l can be derived by 
using Eq. I)A.5(I and the orthonormal condition 



(A.7) 



dt^k{c,t)^i[c,t) = Ski- 

) 

Another important relation is 
2i''RÍ'^{c,l)Sok{c,x) = dye'^^ySok{c,y), 



(A.8) 



from which the Fourier transformation, Eq. (|23|) can be 
derived. From Eqs. (|23|l . and (|A.fi|l . we have 



Xfe(c)$fe(c,í7) 



T/2 



T/2 



dt'^k{c,t)e 



(A.9) 
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